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It is argued that the (traditional) global level statistics which determines localization and coherent 
\Q ' transport properties of disordered systems at zero temperature (e.g. the Anderson model) becomes 

inappropriate when it comes to incoherent transport. We define local level statistics which proves to 
be relevant for finite temperature incoherent transport and optics of one- dimensional systems (e.g. 
molecular aggregates, conjugated polymers, etc.). 
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I. INTRODUCTION 



o 
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Level statistics has extensively been studied in the connection with localization and transport properties of dis- 
i | ordered systems (see e. g. Ref. Q and references therein). In particular, nearest level spacing (NLS) distribution 
function P(S) has a clear physical meaning and its relationship with localization properties is easily established: in 
the localized phase, energy levels are uncorrelated and can be infinitesimally close in energy (the corresponding wave 
+J , functions can be localized far away from each other with a vanishing overlap), which gives rise to the Poisson NLS 
distribution. On the other hand, in the extended phase wave functions overlap well, which leads to strong level 
repulsion and results in the Wigner-Dyson NLS statistics (time-reversal symmetry is assumed hereafter). Various 
quantities associated with P(S) were proposed as scaling variables to analyze metal-insulator transition point, critical 
exponents, mobility edges, and other characteristics relevant for coherent transport (see e. g. Ref. and references 
therein). In the case of incoherent transport a quasi-particle can hop between different localized states because of 
interactions with environment (e.g. with a thermal bath as in the case of the phonon-assisted hopping). The rate (or 
probability) of hopping between two states is then proportional to the square of the matrix element of the interaction 
calculated between wavefunctions of the states. Hopping probability between states with vanishing overlap is therefore 
negligible (localized phase is considered hereafter); scattering rate is relatively large only for the states that overlap 
■ well. The latter introduces correlation between pairs of levels that are relevant for hopping. No such correlation is 
t**»« | accounted for in the traditional global level statistics, we therefore define local level statistics that takes the correlation 
i into account and proves to be adequate for incoherent transport in ID disordered systems. 

In the next section, we introduce a sample model of phonon-assisted hopping, section 3 presents a definition of the 
local NLS statistics in ID, numerical results are discussed in section 4, while section 5 concludes the paper. 



II. DIFFUSION MODEL 



In this section we present an example model of quasi-particle diffusion in a random uncorrelated potential {e„} 

mediated by weak coupling to a phonon system. The quasi-particle spectrum and wavefunctions {E u , are defined 

^ ■ as: 
O 
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H* V = E V * V , H = J2^n\n){n\-t^2(\n){n + l\ + \n + l){n\) (1) 



On-site energies e n are randomly distributed within the box [-W/2, W/2]. All wave functions are localized at W ^ 0; 
incoherent transport can then be modeled by means of the Pauli master equation for the populations P v of the 
quasi-particle states: 

P V = R V - Y V P V + J2( W »» P » ~ W^P V ) , = S(|av|) G(^) I„ M . (2) 

where i?„ is the source term, r„ is the total decay rate of state v (and describes the drain) , W vil is the scattering rate 
from the state to the state \v) 1 which is due to the weak coupling to a phonon system. The rates W V)l are taken to 
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be proportional to the one-phonon spectral density S(u>) and the probability overlap integral I„ M = J2 n *L*ct ( see 
e.g. Refs. 0, H and p| for details), uj^ u = E^ — E v , and G(uj) = n(u) if uj > while G(w) = 1 + n(— w) if oj < 0, with 
n(w) = [exp(w/T) — l]^ 1 being the mean thermal occupation number of a phonon mode with energy oj (the Boltzmann 
constant is set to unity while T is the temperature). The model proved to be relevant for J-aggregates 0, 0,0)0 • 

III. LOCAL LEVEL STATISTICS IN ID 
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FIG. 1: Wavefunctions and a diffusion trajectory: baselines of each wavefunction represent the energy of the state, wave- 
functions are in arbitrary units. Wavefunctions forming the Lifshits tail of the density of states [E v < Eb = —2t, Eb is the 
bare band bottom [at W = 0]) are localized at segments of typical size N* = N*(W). Dashed line — a typical quasi-particle 
diffusion trajectory: because the hopping probability is proportional to the the overlap integral I M „ the particle hops mostly 
over well-overlapping states. AE is the smallest energy separation between two distant states, which is relevant for global level 
statistics (limjv_^oo AE — 0), £j is the energy spacing for the hop i with the typical spacing e* such that limjv^oo e* = const. 



A typical realization of wavefunctions of the Hamiltonian JTJ for a chain of size N and a possible diffusion trajectory 
is shown in Fig. f . Note, that the overlap integral 1^ that enters the hopping rate selects hops between states 
that overlap well (see also the figure caption for comments). The corresponding energy spacing at hop i between such 
correlated levels is a. To account for the correlation we propose the following procedure that is relevant for ID case: 
(i) from the full sorted energy set {E„} select energies within the window of interest (e.g. (— oo, E c ] where the upper 
cut-off energy can reasonably be set as E c ~ Eh or E c ~ Eh + T, Eh being the bare band bottom [at W = 0]), (ii) 
sort the selected energy set in the averaged radius-vector for the state u, x v — ^2 n \^ vn \ 2 n (in ID i is a number 
and such sorting is well defined); let {Ei : Xi+i > Xi] be the x-sorted set, then (iii) construct a new energy sequence 
{£„ : £ n = Y^l s ii s i — \Ei+i — Ei\}, i. e. new energies are cumulative sums of absolute spacings between x-sorted 
energies. The set {£i} takes into account the above-mentioned correlations because sequential states from the set 
{£ n } are always close in r-space, so that overlap integral between them is large. In the next section we calculate the 
local NLS distribution functions Pe (S) for x-sorted energy sequences and discuss the results. 



IV. NUMERICAL RESULTS AND DISCUSSION 



Fig. 2 demonstrates the local NLS statistics results calculated for N = 4096 and averaged over 1000 realizations of 
disorder. Left panel shows the dependence of the mean NLS on the disorder magnitude W . In contrast to the global 
NLS statistics with vanishing mean spacing in the thermodynamic limit, the mean local NLS is independent of the 
system size. Calculated for E c = Eh it determines the characteristic energy scale for the states at the Lifshits tail 
of the density of states, e.g., the diffusion activation energy at low temperature (the latter was demonstrated for the 
case of the Frenkel exciton diffusion, see Ref. Q, where the statistics of the hidden energy structure^ was used; such 
statistics can be viewed as a limiting (less general) case of the proposed local NLS statistics) . The characteristic energy 
e*(W) of the local NLS for E c = Eb is in good agreement with the energy spacing in the hidden energy structure 0- 
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In order to analyze universal correlations between energy levels and compare NLS distribution functions to the 
standard ones (Poisson and Wigner-Dyson) it is necessary to unfold the raw energy sequence {£„} (see Ref. @ for 
details). Right panel of Fig. 2 shows the distribution functions Pe c (S) calculated from the unfolded energy sequence 
for W = 0.3i and different cut-off energies E c . Poisson and Wigner-Dyson distributions are shown for comparison. 
When E c is well below the bare band edge Eb = —2t (E c is deep in the Lifshits tail) the distribution tends to the 
Poisson distribution and demonstrates almost no level repulsion. In this case the states are localized far away from 
each other and there is a vanishing overlap and correlation between them. As the cut-off energy approaches the bare 
band bottom Eb clear sign of level repulsion appears in the local NLS distribution function. However, the distribution 
does not tend to the Wigner-Dyson one upon increase of E c , rather it tends to some limiting distribution. The latter 
is not surprising because although the ai-sorted sequence is correlated it is a sequence of localized states and no true 
Wigner-Dyson statistics can be expected. The reason for this is the following: it turns up that there is a non-negligible 
amount of neighboring segments of the chain with very similar fluctuations of the random potential (averaged over 
the typical localization volume of the wavefunction, N*); energies of the states localized at such fluctuations are 
therefore very close to each other. Such pairs of fluctuations are similar to the system of two identical potential wells 
or a two-atom molecule; as the result, the " eigenstates" of such pairs of fluctuations remind bonded and anti-bonded 
states of two-atom molecule and can be very close in energy (subject to the wavefunction overlap). Positions of both 
states are also very close to each other, so they almost always follow each other in the cc-sorted sequence making even 
the local NLS statistics less repulsive. Quantitative contribution of such pairs of states to the NLS statistics will be 
analyzed elsewhere. 
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FIG. 2: Left panel — the dependence of the characteristic energy scale e* on the disorder magnitude W for E c — Eb. Right 
panel — local NLS distribution function Pe c (S) for the unfolded i-sorted sequence compared to the Poisson distribution Pp(S) 
(completely uncorrelated levels) and the Wigner-Dyson one Pwd{S) (correlated levels). 



The dependence of the local NLS statistics on the cut-off energy allows for definition of the diffusion mobility edge 
in a similar way the conventional mobility edge that separates extended and localized phases is defined. Although no 
scaling variables or rigorous argumentation is available in the present case (to the best of the author's knowledge), 
by analogy, we formally define the diffusion mobility edge as the energy E^ such that: 



PeASo))- 



Pp(S ) + P WD (S ) 



= 



(3) 



where Pp{S) and Pwd{S) are Poisson and Wigner-Dyson distributions respectively, and the spacing So ~ 1. In 
practice, we averaged E^ over various Sq : 0.1 < So < 3. For all values of disorder magnitude the diffusion mobility 
edge coincided with the bare band bottom Eb = —It. In order to interpret the later result it is useful to compare the 
number of states below the mobility edge to the number of the Lifshits tail states N/N* (see Fig. 1), it turns up that 
for all values of W 
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J p (E)dE*—, p(E) = (S28(E-E v 



(4) 



The diffusion mobility edge has therefore a clear physical meaning: it separates the spectral region with states that 
do not overlap with neighbors (hopping rates between such neighbors are vanishing) from the energy region where 
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wavefunctions begin to overlap well: higher (band) states cover one or more tail states and overlap well with them 
and between each other as well (see Fig. 1). In the latter case hopping rates and diffusion build up drastically. The 
importance of two-step hops via higher states was discussed in detail in Ref. Q in connection with the Frenkel exciton 
diffusion. 



V. CONCLUSIONS 



We have defined local NLS statistics for ID disordered systems. The local statistics grasps correlation between levels 
that are involved in hopping transport and is generally more repulsive than the global one. The local NLS distribution 
in the localized phase changes from the Poisson one to the distribution that demonstrates clear level repulsion, 
depending on the spectral region and the energy window. Monitoring the change we determine the diffusion mobility 
edge that has a clear physical meaning: it separates the states with vanishing probability of hopping between neighbors 
(slow incoherent transport) from the states with large inter-state hopping rate (fast transport regime). Using similar 
argumentation local level statistics can be defined for higher-dimensional systems for which the common mobility edge 
can be defined. The mobility edge determined by the local statistics would however be different from the common one 
because the former is always the diffusion mobility edge. Analysis of higher-dimensional cases is beyond the scope of 
this paper and will be discussed elsewhere. 
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